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Abstract 

We extend the formalism of integrable operators a la Its-Izergin-Korepin-Slavnov to matrix- 
valued convolution operators on a semi-infinite interval and to matrix integral operators with 
a kernel of the form E% ^J^ 2 ^ thus proving that their resolvent operators can be expressed 
in terms of solutions of some specific Riemann-Hilbert problems. We also describe some ap- 
plications, mainly to a noncommutative version of Painleve II (recently introduced by Retakh 
and Rubtsov) , a related noncommutative equation of Painleve type. We construct a particular 
family of solutions of the noncommutative Painleve II that are pole-free (for real values of the 
variables) and hence analogous to the Hastings-McLeod solution of (commutative) Painleve 
II. Such a solution plays the same role as its commutative counterpart relative to the Tracy- 
Widom theorem, but for the computation of the Fredholm determinant of a matrix version of 
the Airy kernel. 
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1 Introduction and results 

The paper aims at extending the general theory of integrable operators of ItsTzergin-Korepin- 
Slavnov (IIKS for short) [11] to operators of "Hankel" form (see below). Leaving aside for the time 
being any analytical consideration, the issue is the study of integral operators on L 2 ( r y + ,C r ) with 
a kernel of the following form 

MA,/i) = — , A,ptG7 + (1.1) 

A ■+■ fX 

where 7+ is a contour contained in a half-plane of C (so that the denominator does not vanish) 4 and 
the matrices Ej : 7+ — > Mat(p x r, C) are suitable (analytic) functions. These operators are related 
via a Fourier transform to (matrix) convolution operators on R+ as pointed out in Section 2: our 
primary focus shall be the construction of a suitable Riemann-Hilbert problem for computing the 
resolvent operator S = —K, o (1 + K,) -1 . The knowledge of the resolvent operator allows to write 
variational formula; for the Fredholm determinant of the operator Id + /C : L 2 (j + , C) — s- L 2 (-y + , C r ) 
via the well-known variational formula 

91ndct(Id + /C) =Tr((Id + 5)o9/C) . (1.2) 

4 The contour could be -for example- R+ provided that the matrices in the numerator yield _E^(0)i?2(0) = 0. 
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The situation is closely related to the IIKS theory mentioned above (with the tensorial extension 
explained in [8]), which we briefly recall: let £ C C be a collection of (smooth) contours and let 
f,g:E-) Mat(g x n, C) be smooth (analytic) functions on 7 subject to the condition 

f T (A)g(A) = 0, Ae£. (1.3) 

Consider the integral operator N : £ 2 (£, C n ) with kernel given by 

_ (1 .4, 

A — n 

Then the resolvent operator 1Z = N o (Id — TV) -1 has a kernel (denoted with the same symbol 1Z) 
of the form 5 

KKv) = PW T M e - T MM> (1.5) 
A — (i 

where 0(A) is the q x q matrix bounded solution of the following Riemann-Hilbert problem 

0(A)+= 8(A)- (1, - 2*7rf (A)g T (A)) 
9(A)= 1, + (D(\- 1 ) , A ^00 (1.6) 

Furthermore the solution of the RHP (1.6) exists if and only if the Fredholm determinant det(Id — N) 
is not zero. 

The operator (1.1) is not immediately of the form (1.4) and hence the IIKS theory is not directly 
applicable. Nevertheless the former situation is amenable -not surprisingly- to the latter (see also 
[23]). In fact one could observe, for example, that K? is an operator of the form (1.4) 

K 2 (A,,,)= jl— (Ef - ff 2 (A)J5 2 M) 

BM luM>M®*M, (1J) 

This observation shows that the IIKS theory is relevant also to the study of operators of the form 
(1.1): however it is not practical to use (1.7) as a starting point for the analysis as this route is 
impervious and is not the one we follow. We provide a direct treatment of IC as well as JC 2 in a 
unified fashion; the RHPs that are relevant are specified in Problems 3.1, 3.2 (please refer to the 
statements there) for two matrix functions T, 3 of size 2r x 2r. The two problems are intimately 
related to each other in that the jump conditions are identical while only the asymptotic behavior 
at A = 00 for r, S differs. The solubility of the Riemann-Hilbert Problems 3.1, 3.2 is equivalent 



3 The superscript to a matrix denotes the inverse transposed matrix. 
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to the non vanishing of the Fredholm determinants of the operators Id 7 — JC (Thm. 3.1) and 
Id 7+ + JC (Thm. 3.2), respectively, which follows from IIKS theory; we thus obtain the formula for 
the resolvents of JC and JC 2 (Theorems 3.2 and 3.1) 





2fi [El(\),O pXr ] T T (A)T~ T (/i) 
- JC o (Id 7+ + JC)-i(\, fx) = 



rxp 
E2^) 



/C 2 o(Id 7+ -K?y\\ii) = [Ef(X),O pXr 



5 r (A)5-^) 
A — fi 



O r xp 



(1.8) 



(1.9) 



The knowledge of the resolvent operator allows to write variational formula? for the respective 
Fredholm determinants: however one may bypass formula (1.2) and write the variational formula; 
directly in terms of the solution of the respective RHPs (Thm. 4.1 4.2) using the ideas in [2] 

<91ndct(Id 7+ -/C 2 )= f Tr^S'.dMM- 1 ) — (1.10) 

1 r a \ 

<91ndct(Id 7+ +JC)= - I Tr (FZ^dM M' 1 ) — (1.11) 

M(A) := l 2r - 2inE 1 (X)E 2 (\) T ( X , + « cr + + ® <r_ ) , 7 - := -7+ (1.12) 

where / is derivative w.r.t. A, 9 denotes any variation of the symbols Ej and <r+ ( cr -) denotes the 
2x2 matrix with just one non-zero entry on the upper right corner (lower left corner). 

In the second part of the paper we provide some applications to the study of matrix convolution 
operators; our example of choice is a matrix version of the (scalar) convolution operator by the Airy 
function [6] 

Ai s :L 2 (R+) -> L 2 (R+) 

f(y) ^ (Ai s f)(x) := / Ai(.T + y + 2s)f(y)dy (1.13) 



The Fredholm determinant of the operator Id — Ai s /2 is known to yield the Tracy- Widom gap 
distribution -Fi(s) for the GOE [6] and -on the other hand- the Fredholm determinant of Id — Ai 2 
yields the distribution F 2 (s) for the GUE [17] ; in fact it is well known that the kernel of the square 
of the Airy-convolution operator is the celebrated Airy kernel 

M 2 (x 7 y):= [ Ai(x + z)M(y + z)dz = Ai ^ Ai '^ ~ A[ (yW(^) = , y ( } (1 u) 
Jr + x-y 

F2 is expressed in terms of the Hastings-McLeod solution [9] to the second Painleve equation [22] 
while F\ can be expressed in terms of the Miura transform of the same transcendent. Alternatively 
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(and equivalently) -Fi(s) can be expressed in terms of the unique solution of the Painleve XXXIV 
equation with a certain prescribed asymptotics 6 . 

*«- (- jy - .«,)■*) , ^ „m = i„ « w ,i, 5 , 

iq(s)=exp(-i / u(x)dx) (f 2 ( s )V (1.16) 



2„ 

u"(s)= 2u(s) 3 + su(s), u(s)~Ai(s), s — > +00. (1.18) 
w"'(s)=12w(s)w'(s)+2w(s) + sw'(s), w(s) ^Ai'(s), s — ► +00. (1.19) 

w(s)= \u\s) - \u'{s) (1.20) 

where (1-20) is the usual Miura transformation between solutions of modified KdV and KdV equa- 
tions. Moreover we refer to (1.19) as the PXXXIV equation since, up to rescaling, this is the same 
as the derivative of equation (30) in [4] (see also [16]). 

The noncommutative analog of the whole preceding discussion arises in the study of a ma- 
trix version of the Airy-convolution operator (see Section 5) which we have picked as exemplary 
application: 

{Aisf)(x):= f Ai(x + y;s)f(y)dy (1.21) 
Ai(x; s):= [c jk Ai(x + 8j + s k )] j k . (1.22) 

Here the matrix C = [cjk]j k is an arbitrary r x r matrix with complex entries (in general) and the 
dependence of Aig on C is considered as parametric (and it is understood in the notation). The 
kernel of the square of this matrix-kernel docs define a probabilistic model because it is a totally 
positive kernel on the configuration space {1, . . . ,r} x K as shown in Thm. 5.2. The analysis of 
Aig and Aig is then related to certain noncommutative analogs of the aforementioned Painleve 
equations and particular solutions thereof; in particular the Fredholm determinant of Ai\ is related 
to the noncommutative (matrix) Painleve II equation 7 

r o 

D 2 C/(s) =4(slI(s) + U(s)s) + 8U 3 (s) , s := diag (si, . . . , s r ) , D := ■ ( L23 ) 

.7 = 1 Sj 



6 The uniqueness of the solution w with the prescribed asymptotics is easily deduced from the uniqueness of the 
Hasting-Mc Leod solution u of PII. 

^Equation (1.23) reduces to (1.18) in the scalar case r = 1 with the change of variable x = 2s. Also, the matrix 
U(s) in the body of the paper shall be denoted by 0i(s). 
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which appeared recently in [19]: that paper provided special solutions in terms of quasideterminants 
[7] in a more general context of noncommutative rings, but not a Lax-pair representation or a 
connection to Ricmann-Hilbcrt problems or Frcdholm determinants. The isomonodromic approach 
to the above equation yields a Lax pair representation contained in Section 5.1 and particularly 
Lemma 5.1. 

Of greater interest is the fact that the particular solution that is involved in the computation of 
the Fredholm determinant of Ai\ enjoys the same smoothness properties for s £ R r as the Hastings- 
McLcod solution. More precisely we prove (Prop. 5.1) that there is a unique solution of (1-23) with 
the asymptotic 

U M (s) = c M ki{s k + s e ) + (■/Se-^ {2S - 2m) ^j , S := ^ ]T Sj , m = max \ Sj -S\, S^+oo . 

Additionally, this solution is pole free for s € R r if the maximal singular value of the matrix 
C = [cki] is one or less 8 , a condition which is sufficient if C is an arbitrary complex matrix and 
becomes also necessary if C is Hermitean (Thm. 5.1 and Thm. 5.3). The analog of the third-order 
ODE for Fi is now a system with noncommutative symbols (Thm. 5.4) that can be reduced to a 
fourth order matrix ODE (Remark 5.6) and only in the scalar case is further reduced to an ODE 
of the third order. The Fredholm determinant of Aig (the analog of Fi) is then computed in terms 
of the relevant solution in Corollaries 5.2, 5.3. 

2 Matrix convolution operators on a semi— infinite interval. 

Given a function C : R — > Mat(r x r) decaying sufficiently fast at infinity let's consider the 
convolution operator C acting on L 2 (R + , C) as follows: 

/>oo 

(Ccp)(x) = C(x + yMy)dyeL 2 ([0 1 ^) 1 C r ) (2.1) 
Jo 

Our aim is to study the Fredholm determinants det(Id + C) and det(Id — C 2 ) 9 . Here C (and hence 
the determinant) may depend on some parameters not explicitly indicated here (see below). Such 
type of determinants appears in many applications; just to cite two of them let's recall the Dyson 
formula [5] in the inverse scattering for the Schrodinger operator and in the integral formula of the 
Tracy- Widom distribution for GOE found by Ferrari and Spohn [6] (see below) for r = 1. 

8 The singular values of a matrix C are the (positive) squareroots of the eigenvalues of C^C: they coincide with 
the absolute values of the eigenvalues of C if it is Hermitean (or more generally normal). 

9 Of course the sign in the expression det(Id + C s ) is inessential since we can always change C with — C. 
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Remark 2.1 In the inverse scattering theory of the Schrddinger operator and other applications 
the Fredholm determinant is written as the restriction to [s, oo) of the convolution by C: 

/oo 
C(x + y)f(y)dyeL 2 ([s,^)). (2.2) 

This is identical to the setting above up to translation. In fact it is just enough to redefine C(x) i — > 
C s (x) := C(x + 2s) and make it act on L 2 ([0, oo). 

We will consider functions C(z) that admit the following representation (the factor of —i being 
purely for later convenience) 



C(z) = -i / e^r(/i)d/i (2.3) 

where 7+ stands for a finite union of oriented contours in the upper-half plane with positive distance 
from R and r(/z) is a bounded L 1 ( r y + , Mat(r x r)) function on 7+ (with respect to the arc-length 
measure). This assumption guarantees that C(z) is rapidly decaying at z = +00 G M with a simple 
estimate 10 

\C(z)\ < c - zdist (T+' R ) [ |r(M)||d/i|. (2.4) 

An interesting example is as follows 
Example 2.1 Lei C(z) = — Ai(z) and r = 1; then 

C(z) = -Ai(z + a/2) = ~ [ c l T+»(-+-/2)M d/i (2 .5) 

where 7_|_ is a contour extending to infinity along the directions arg(/i) = ^ ± ^ . TTus example is 
relevant for applications since, as we have written in the introduction, the Fredholm determinant of 
the corresponding convolution operator is equal to the Tracy- Widom distribution for GOE, namely 

Fi(s) = det(Id + C). 

We would like to transfer the study of the Fredholm determinant of C on L 2 (M + ) to the study 
of a Fredholm determinant of an operator in L 2 (7 + ); this is accomplished hereafter. 



10 The symbol |r| on a matrix stands for any norm on the matrices, for example the Hilbert-Schmidt norm or the 
supremum of the absolute values of the entries. This is so not to overload the notation when considering norms in 
some L p . 
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Proposition 2.1 Let C(z) as above, with r(/i) = E\(ij)E% (p) and Ej eL 2 fl L°°(^ + , Mat(r x p)); 
then the operator C is of trace-class on L 2 (R + ,C r ) and also 

det(Id L 2 (R+) + C) = det(Id^2 + K) (2.6) 

where Ti 2 . is the Hardy space % 2 ® C r (i.e. the unitary image of the Fourier-Plancherel transform 
of L 2 (R+, C r ),) and where JC is the integral operator on H 2 C L 2 (R, C) with kernel 

Proof of Prop. 2.1. By Paley- Wiener theorem, L 2 (R + , C r ) is unitarily equivalent under Fourier- 
Plancherel transform T to the subspace T-L 2 := T-L 2 (g>C r , with Ti 2 the Hardy space of the upper half 
plane. Hence the convolution operator acts as follows 

/>oo />oo /> 

1>(x) := (Cip)(x) = C(x + y)(p(y)6y = -i dy d£e^+^r(£My) = 

JO Jo J~/ + 



= -iV2n /d£e l3 *r(£)(7»(0 

so that (note that the x-integral below is convergent because £ G 7+ C C+) 

(TV-) (A) = -= / e^^dx = -i / dxe iA * / d£e**r(0(7»(0 = 
V 27T Jo Jo j 7+ 

-/^x^^o-/*^"'^- <2 - 8) 

We note that for a function in H 2 like /(//) := T<p(fi), the evaluation at a point £ G C+ can be 
written as 

m = I m . r 1 . ( 2 - 9 ) 

which is Cauchy's theorem. We shall thus define 

Z T := T~ X CT (2.10) 
(the reason for the transposition is solely for later convenience) with kernel given by 

^/(A) ^^^v - £/(A) = / d, / di (2.H) 

Finally, since the Fourier Plancherel transform from L 2 (R + , C r ) to "H 2 is an isometry, the respective 
Fredholm determinants are equal (if they exist). We note that tC extends to an integral operator on 
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the whole of L 2 (R, C r ) with the same kernel: this extension automatically annihilates the orthogonal 
complement of H 2 in L 2 (R, C), which is seen by closing the ^-integral with a half circle in the 
lower half plane and then invoking Cauchy's theorem. We will understand this extension in what 
follows. 

Therefore to conclude we need to show that C and K. are trace-class. By the unitary equivalence 
given by the Fourier-Plancherel transform it suffices to show that JC is trace class: we shall present 
K, as the product of two Hilbert-Schmidt operators, thus proving it of trace class. To this end recall 
that r(fj) = £'i(/i)£ , J , (/i); thus 

KfW = /dM / d£ ~r r , —z-s =C 2 o Cif(X) . (2.12) 

where the two operators are defined as follows 

Cr. %l c L 2 (R, C r ) -> L 2 ( 7+ , C p ) 

/^(d/)(€):=^(0 / /t%t (2-13) 
Jk 2itt(^ + £) 

C 2 :L 2 ( 7+ ,C7)^H 2 

^(« )M:= ^^M (2 , 4) 

We embed % 2 and L 2 ( 7+ , C p ) as subspaces of L 2 {R U 7+ , C r+P ) as 



/(A)X K (A) 



0? 



L 2 ( 7+ ,C P ) 9 /ih-> 



(2.15) 



(they are then orthogonal but not complementary) and think of Cj as extended to the whole 
L 2 (R U 7 +,C' +p ) in the trivial way (i.e. acting like zero on the orthogonal complements of the 
H 2 , L 2 ( 7 +,C P ), respectively). Analogously we extend trivially the action of K to this enlarged 
Hilbert space. Then it is promptly seen that they both are Hilbert Schmidt in L 2 (lU 7 +,C +p ) 
because 

w <+ °° <2j6) 

thanks to our assumption that (the entries of) Ej are all in L 2 ( 7+ ,C). Thus K, : % 2 — > "H 2 is trace 
class, so is C and their determinants are the same. Q.E.D 

Recall that if A : Hi —> H.2 and B : H2 — > "Hi are (bounded) operators between Hilbert spaces 
and both AB, BA are trace class then (see for instance [20]) 

Tr(Bo4) = Tr(4oB). (2.17) 
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Composing the operators Cj in the opposite order we obtain an operator on L 2 (7 + ,C P ) as follows 

The ^-integral can be closed with a big circle in the upper-half plane, thus picking up the residue 
at £ = A to give 

(Ci o C 2 /)( M ) = El^) I d\ E2{X)fiX) =: (£/)(/,). (2.19) 

J 1+ A -\- fJi 

Renaming the variables we obtain that 

^.smsd. (2 .20) 

A + H 

The operator JC is also trace class because the composition of two Hilbert-Schmidt operators in 
L 2 (7 + UK, C r+P ); clearly it defines an operator on L 2 (-f + , C p ) since it acts trivially on its orthogonal 
complement (by construction). In particular 

Corollary 2.1 The Fredholm determinants of C : L 2 (R + ,C r ), K. : U 2 r -> H 2 r and K, : L 2 (j+,C P ) -> 
L (7_|_,C P ) are all equal. 

Given the fact that the operator JC (with kernel (2.7)) and JC (with kernel (2.20)) have the same 
Fredholm determinant, we shall continue our discussion by focusing on the latter. In the sequel 
we will simply analyze operators with kernels as in (2.20) and forget about their origin as Fourier 
transform of convolution operators. 

Remark 2.2 It is worth mentioning that operators with kernel (2.20) withy = r = 1 (and Ei = E2) 
belong to the same class considered in [23]. Using our formalism it is possible to re-derive the 
connection between the Fredholm determinants of these operators and the mKdV/KdV hierarchies. 

3 Riemann-Hilbert problems with different asymptotics and 
their mutual relationship 

Given a kernel JC(X, fi) as in (2.20) corresponding to the operator (denoted by the same symbol) 

/C:L 2 ( 7+ ,C p )^i 2 ( 7+ ,C p ), JC(X^) (3.1) 

A + (j, 



we construct two related Riemann-Hilbert problem on the collection of contours 7 
7_ := —7+) and with jump matrix 



= 7+ U 7_ (here 



M(X\- lr -2^r(A) X 

r(A) = Bi(A)E| , (A) e Mat(r x r) , r(A) := r(-A). 



(3.2) 
(3.3) 



where \ x denotes the indicator function of the set X. Here and below we denote with erj, i = 1, 2, 3 
the Pauli matrices 





' 


1 " 




i ' 




' 1 


(7i := 


1 





cr 2 := 


-i 


0-3 := 


-1 



cr + := (5ii5j2)i j=ifi and er_ its transpose. Furthermore we shall set 

(7 fe = l r ®cr k , k = ±, 1,2,3, 



(3.4) 



where by the tensor notation can be taken to mean the matrix of size 2r x 2r split into 2x2 blocks 
of size r x r. 

Note that the jump matrices M(A) on 7 := 7+ U 7_ satisfy 



M(-A) =ctiM(A)cti. 



(3.5) 



We are going to formulate two Riemann-Hilbert Problems (Problems 3.1, 3.2) and we will show 
how they are related between themselves (Prop. 3.2) and how they relate respectively to the two 
Fredholm determinants det(Id 7+ + K.) and det(Id 7+ — K?) and the resolvents of the respective 
operators. In the sequel we shall assume that Ej(X) are smooth (beside the already imposed 
conditions Ej G L 2 (j + ) n L°°(7+)®Mat(r xp)). 



Problem 3.1 Find the sectionally analytic function H(A) G GL(2r,C) onC \ (7+ U 7-) such that 
(with M(A) given in (3.2)) 

5+ (A) = 5_(A)M(A) Ae 7+ U 7 - (3.6) 

5(A) = lar + ^ + + • • • , A -> 00 . (3.7) 
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Problem 3.2 Find the sectionally analytic function T(A) G GL(2r, C) on C \ (7+ U 7_ 
(urai/i M(A) given in (3.2)) 


) such that 


r+(A) 


= r_(A)M(A) 


A e 7+ U7_, 


(3.8) 


T(A) 


= L(A) ^l 2r 4 




(3.9) 


rCAjL-^A) 

r(-A) 


= 0(1) , A - 

= r(A)Si 

= ai»tT3 


-> 


(3.10) 
(3.11) 
(3.12) 


where the matrix L(A) is defined as 


follows 






L(A) := l r »L(A) = 


—iXl r i\l r 


> £ W - [ 4a I 


(3.13) 



The validity of the asymptotic expansions near infinity needs additional conditions on the jump 
matrices if some component of 7+ extends to infinity (if this happens we assume that these com- 
ponents extend to infinity along asymptotic directions) . A sufficient condition, which we hereby 
tacitly assume, is that r(A) = 0(|A|~°°) as |A| — > 00 along any such component and extends to an 
analytic function on an open sector containing the direction of approach in such a way that the 
same asymptotic holds. 

It is clear that the two Problems 3.1, 3.2 are closely related and the remainder of this section 
is devoted to explaining their mutual relationship. It is a straightforward result that, if a solution 
of Problem 3.1 exists, then it is unique. The uniqueness of the solution for the Problem 3.2 comes 
from the following proposition. 

Proposition 3.1 Let T(A) be a sectionally analytic function that solves the RHP (3.8, 3.9, 3.10, 
3.11). 

If a solution exists then 

1. detr(A) = {2i\) r ; 

2. Any matrix T = (l 2r + c®G-)T solves the same RHP with c € Mat(r x r), constant; 
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3. Any solution has an expansion where the terms Tj in (3.9) have the symmetry 



r(A) = L(A) [Ur + Y. 



X k 



) 



r, = (-1)^1^ 



(3.14) 



and hence 



I A., = a 2j »l2 + ''j, T| , T 2 j+i = "2j ■ i + b 2] +i»o-2 , € MiU ; r x r) (3.15) 

4- If we additionally require the condition (3.12) T\ = ai®(73 (for some constant matrix a\) then 
the solution is unique. This solution will be referred to as the gauge-fixed solution. 

Proof of Prop. 3.1. 1. It is clear that the determinant has no jumps because the jump matrices are 
unimodular. Moreover, from detL = (2iX) r (see (3.13)) and (3.9) we have detr(A) = (2i\) r (l + 
©(A" 1 )). Finally from (3.10) we have detT = O(dctL) = 0(X r ), X -> 0, and hence it must be 
dctr = (2iX) r . 
2. We note that 



and hence the multiplication on the left by a constant matrix of such a form does not change the 
form of the asymptotic expansion and does not change the jump conditions. This proves the second 
point. 

3. The statement is obvious once one notices that L(— A) = L(A)(?i. 

4. Suppose r is another solution satisfying the same requirements and denote by aj,bj the 
coefficients in its expansion as per (3.15). By point 1, any two solutions have the same determinant; 
the ratio 



must be a holomorphic matrix function on C \ {0}. However, from the condition (3.10) we see that 
actually S(X) must be analytic at as well. Looking at the behaviour at infinity of T and T one 
hnds by a direct computation that S(X) is bounded and 




(3.16) 



5(A) :=r(A)r-!(A) 



(3.17) 



5(A) 



lr 

ia\ — ia\ l r 



(3.18) 



Suppose now that a\ ^= ai; then 



r(A) 



l r 

ia\ — ia\ l r 



r(A) 



(3.19) 
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But then one sees by direct matrix multiplication that 61 should equal ai 2 ai which violates the 
normalization b\ = 0. This proves uniqueness. Q.E.D 

In Proposition 3.2 and Proposition 3.3 we study the relationship between the Riemann-Hilbcrt 
problems 3.1 and 3.2: in particular we shall see that they are not equivalent, in the sense that 
if Problem 3.1 admits a solution then so does Problem 3.2 but, in general, not viceversa. We start 
by observing that the symmetry (3.5) for the jump matrices implies the same symmetry for 3 

3(-A) = ?iS(A)3=i (3.20) 

which in turns implies the following form for the coefficient 3j in (3.7) 
00 „ 

3(A) = l2r + ^ Tx: i s 2j+i = a 2 j+ie:(J3 + /32j+i®02 , 3 2 j = a 2 j®l 2 + /32j® o"i . (3.21) 
fe=i 



Proposition 3.2 Let 3 be the solution of Problem 3.1; then the solution of Problem 3.2 is 



r(A) = 



-iXl r - 20! iXl T - 2/3i 



3(A) 



(3.22) 



iwii/i /?i as m (3.21). In addition the coefficients of the expansions for T (3.15) and 3 (3.21) satisfy 



01 = ai - ipi , 
fl2j = «2j - («2j-i - ifaj-i) , 



&i = (3.23) 

b 2 j = lh 3 - ifii (asy_i - »Ay-0 (3.24) 



Proof of Prop. 3.2. Since T and 3 have the same jumps we must have T(A) = i?(A)S(A) for some 
i?(A) at most polynomial. From the symmetries we must have R(—X) = R(X)ai and dot R = (2iX) r . 
The expansion of T and 3 at infinity forces R to be of the form 



R(X) = 



— «Al r + 2ic iXl r + 2ic 



On the other hand, as we presently show, the gauge fixing (3.12) determines C; indeed 



L 1 R — I2,' 



A 



-1 -1 
1 1 



1 . 

= J-2r + T (— + ic»<7 2 ) 
A 



and therefore in the expansions of T and 3 and matrix multiplications we have 

a2j+l®<73 + &2j+l®CT2 V"^ a 2j®l-2 + &2j®CTl 



J=0 



1 5^ 

J=l 

(«i — c)«CT 3 + (/?i + ic)®(T 2 
A + 



L2r 



(3.25) 



(3.26) 



(3.27) 
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(ct2j - c(a 2 j-i - ^2j-i))«'l2 + (Jhj - c(a 2 j-i - i/3 2 j-i)) »<Ti 
. ("2^+1 - ca 2j + c/3 2 j) ®tr 3 + - ic^ 2 j + «ca 2 j) ®er 2 , , 

+2^ ^ ( } 

The gauge fixing (3.12) mandates 61 = (i.e. the coefficient matrix of a 2 in the term A -1 must be 
absent) so that we must have c = ifii and equating the coefficients of the expansion above implies 
(3.24). It only remains to show that i?(A)S(A)L _1 (A) is bounded at A = (condition (3.10)). Since 
L,- 1 = l^L- 1 and £ _1 (A) has only simple pole at A = 0, then i?(A)S(A)L _1 (A) = f + 0(1) as 
A — >• 0. On the other hand the symmetries imply i?(A)S(A)L~ 1 (A) = R(— A)H(— A)L _1 (-A) and 
hence Co = — Co so that Co = 0. Q.E.D 

Proposition 3.3 Let T(X) G GL(2r, C) fee i/ie solution of Problem 3.2 and denote the r x r blocks 
of r fey Py, z, j = 1, 2; £/ien £/ie solution of Problem 3.1 for S ea;iste i/ and ora/y ?/ 

detru(O) 7^ . 

Moreover 

01 = -ir n 1 (0)r 2 i(0) = -i lim AS 12 (A) 

Z A— >oo 

Proof of Prop. 3.3. Let T(A) be the solution of Problem 3.2. In particular T(A) is bounded 
everywhere (by definition) and we want now to find a matrix R(X) of the form (3.25) such that 

5(A) := i? _1 (A)r(A) (3.29) 

solves Problem 3.1. It is clear that the jumps will be automatically satisfied and so the asymptotic 
behaviour at infinity. The value of the constant matrix c must be determined by the requirement 
that S is bounded at A = 0. From the symmetry (3.11) we have the matrix equations 

rn(o) = ria(o) , r 2 i(o) = r 22 (o) . (3.30) 

A direct computation yields 

i?~ 1 (A)r(A) = 0(1) , A^0 (3.31) 

T(A) = 0(1) (3.32) 





c 
A 


1 T - 

2iX 


lr 

- 2 


c 
A 


l r 
2i\ . 


) = o, 


=> 


i 

C= ~2 



lr 21 (o) = o,=* c = --rr 1 i (o)r 2 i(o) (3.33) 
c r 12 (o) + ^r 22 (o) = , (3.34) 
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The two equations are the same due to (3.30). Now, if det rn(0) 7^ then the solution for c is as in 
(3.33) and the sufficiency is proved. As for the necessity, if detr n (0) = 0, then the equation (3.33) 
may still be compatible. However this would mean that there arc infinitely many c that solve the 
matrix equation (3.33), which would violate the uniqueness of the RHP 3.1. Q.E.D 

We conclude the section with the following two theorems, which we state side-by-side for the 
sake of easy comparison. 



Theorem 3.1 Let /C(A,/i) be the integral operator on L 2 (j + ,<C P ) with kernel 

A + fi 

Then the resolvent operator — K? o (Id 7+ — /C 2 )" 1 0/ Id 7+ — K? on L 2 ('j + ,C P ) has kernel 

^++(A,Ai) given by 



TZ ++ (X,n) = [£?(A),0 



,h t (a)h- t (m) 



px r J 



Or xp 



(3.36) 



X — fi 

where S is the solution of Problem 3.1 with the jump matrix (3.2). 

If r(A) := i?i(A)£ , J(A) is symmetric, r = r T (for example if E\ = E 2 = E) then the resolvent 
can be written more symmetrically as 

] E T (X)a 2 E(fi) 



n ++ (X,fx)=i[El(X),O p> 



X — fi 



0, 



(3.37) 



The solution to Problem 3.1 exists if and only if the operator Id 7+ — IC 2 is invertible. 



Theorem 3.2 Let K,{X,fi) be the integral same operator as in Thm. 3.1. Then the resolvent 
operator S — —K o (IcL, + /C) _1 has kernel S(X, fi) given by 



2fi [El(X),O pxr ] T T (X)T- T (fi) 



S(X,fi) 



E 2 (fi) 



(3.38) 



A 2 — fi 2 

where T is the solution of Problem 3.2 with the jump matrix (3.2). If r(A) := i?i(A)£'J(A) is sym- 
metric, r = r T (for example if E\ — E 2 — E) then the resolvent can be written more symmetrically 



15 



as 



S(X,li) 



[Ef(X),O pXr ] r T (A)a 2 r(/x) 



Or 



Jrxp 



A 2 



(3.39) 



This solution to Problem 3.2 exists if and only if the operator Id 7+ + K, is invertible 



Proof of Thm. 3.1. We start observing that the jump M(A) in Problem 3.1 can be written 



as 



Af(A)= 1 - 2i7rf(A)g T (A) 
f(A)= 

g(A) = 



" Si (A) " 

O r xp 


x 7+ (A) + 


Oj-xp 

£ X (A) 


O r xp 

_ E 2 (\) _ 




' S 2 (A) " 

O r xp 



X ,_(A), Sj(A) := Ej(—X). 



(3.40) 
(3.41) 

(3.42) 



By the IIKS theory, this RHP is associated to the kernel N acting on L 2 (j + U 7_, C p ) with kernel 
given by 

f T (A)g(/i) _ BfWliW^ Wx 7 _ W + ^(A)£j Wx,. (A)X T+ W 



N(X,n) = 



(3.43) 



X — fx X — fi 

According to the split L 2 (7+ U 7_ ) = L 2 (7+ ) © L 2 (7. ) , using the naturally related matrix notation, 
we can write N as 



N 






G ' 


J" 






£:L 2 ( 7 _,C7)^L 2 (7+ 



^:L 2 ( 7+ ,n^L 2 (7_,CP) (3.44) 



where the operators J- and Q are integral operators with kernels 



EUX)E 2 (Q x ^(X)x^(0 



ET(0E 2 (^_(0x y M 



£-0 



We observe that the kernel of the composition reads 

m)ET(o 



(gor)(x, f i) = El(x) 

Ef(X) 



_ (a-occ-m) 

jg 2 (fliff(fldA 

+ (A+e)(e + ^)i 



^ ^ 2 (m) 



E 2 (/z) =/C 2 (A, m ) 



(3.45) 

(3.46) 
(3.47) 



and hence our task of computing the resolvent of Id 7+ — K? is the same as computing the resolvent 



of Id 7+ — Q o T . To this end we write first the resolvent of Id 



7+U7- 



N (using [11]): 



K(X, fi) = 



f r (A)5 r (A)5-^)g(^) 
A — \i 



(3.48) 



16 



according to the projections in L 2 (j±,C p ): 



[Ej(X),Op X r 

-[Ei (X),0 pxr 



5 r (A)5-^M 

A — fj, 
S T (A)5^( M ) 

X — 



Orxp 



X T+ (A)X 7 >) + [Opxr^iX)} 



Z T (X)E- T (») 



X T+ (A) X ,» + [0 P xr,S 1 T (A)] 



A — /i 
S T (A)S-^M 

X — 



O r xp 



x 7 _(a)x 7 _0) 

(3.49) 



where the four addenda appears in the matrix notation induced by the splitting L 2 ( 7+ U 7_ ) 

L 2 ( 7+ )©L 2 ( 7 _); 



(Id 7+U7 _ 
On the other hand we have 



K)~ 



Id, 



■Ha 



K- 



IcL 



(Id 7+ -QoF)- 



Fo (Id 7 -go J")- 1 Id 



Id 7+ — Q oT 



Id, 



Id, 



-J" Id, 



Id, 



Id, 



Id, 



11— 



so that the entry (1,1) of the equation above gives 

(id 7+ -go ry 1 -- 



Id 



7+ 



Id 7 _ + TZ— 



id. 



Id, 



and the equation (3.49) gives the precise form of the kernel 7£++(A, //). 

In case of symmetry r = r T this form simplifies because S _1 (A) = a 2 S T (X)a 2 (which is proved 
along the same lines as in Thm. 3.2). The statement about the existence is a direct application of 
IIKS theory. Q.E.D 

Proof of Thm. 3.2. The idea of the proof is to reduce as much as possible the theorem to the 
theory of intcgrable operators of Its-Izcrgin-Korepin-Slavnov (IIKS) . We can write the operator as 



lC(X,n) 



E 1 (X) T E 2 ([i) _ (\- f i)E 1 (X) T E 2 (p) 



(3.50) 



A + fx A 2 - fi 2 

We now introduce the coordinate z := A 2 and w := /i 2 . Since 7+ is in the upper half-plane, its image 
under the square map is well-defined and lies in C \ M+ . Since the arc- length of 7+ differs in the 
z-plane and A-plane, we must introduce the square-roots of the Jacobians. The integral operator 
(3.50) reads 



IC(z, w) 



(^-^)E 1 (^) T E 2 (^) 
2(zw)i (z — w) 



E 2 (V^) 



2(z-w) 



(1) 



'3.51) 
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We have to construct the resolvent of Id 7+ + K, = Id 7+ — (— K). We have now an integrable kernel 
in the sense of Its-Izergin-Korepin-Slavnov where the matrices f , g can be chosen as 



(-fC)(z,w) 



f T (z)g(w) 



z — w 



f(z) 



1 



-iy/zEl(y/z) 



We immediately observe that 



f(z) 







rxp 
rX p 



L(^)f (z) , 
= L- T (z)g (z). 



(3.52) 

(3.53) 
(3.54) 



The construction of the resolvent is then associated, in the standard way [11], to the following 
Riemann-Hilbert Problem 



9(z) + = 6(z)_ (l 2r - 2*7Tf (z)g T (z)) , ZG 7+ 

6(z) = l 2r + £>(z~ 1 ) , z^oo 



We can rewrite the jump matrix as follows 



l 2r - 2z^L(VI)(r(VI)«(T+)L- 1 (^) 
r(Vi) :=£?i(Vi)^(>/i) 



(3.55) 
(3.56) 



(3.57) 
(3.58) 



Consequently wc introduce the new matrix 0(z)L(y / z), where li(y/z) := l r ®L(i/z). In order to 
connect with Problem 3.2 we define 



r(A) = e(A 2 )L(A) 



(3.59) 



and we see immediately that f (-A) = f (A)cti. Furthermore f (A)L _1 (A) = 6(A 2 ) = 0(1) as A -» 0. 
Thus T solves Problem 3.2 except for the gauge-fixing (3.12), which we now take into consideration: 
if we denote by a± the (1, 2) block of size rxrin 



e(z) = i 



i 



-<Zi 



+ 0(z~ 2 ) 



(3.60) 



then one verifies by matrix multiplication that the relation between T and the gauge-fixed T is 

(lar + ai»a+)T(z) =T(z) , (3.61) 
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The resolvent operator, according to the general theory, is 



5(A ^ ) = fW(^ w ) g H^ ; 

z — w 

z — w 

f -(A)r-(A)r-(,) g0 (M)^ : 

rX p 



2\A/i 



A 2 - /i 2 

2 M [£f(A),o pxr ] r T (A)r- T (/i) 



A 2 - A j 2 



Finally, note that 



and if r = r T then 



L (A) = —d 2 h(X)d 2 



r t (m) = 77— s : 2 r(^)(T2 



2i/i 



(3.62) 
(3.63) 
(3.64) 

(3.65) 
(3.66) 
(3.67) 



which can be checked by verifying that li^B^Y T (m)^2 solves the same Problem 3.2 and hence 
equals r. In this case the formula for the resolvent takes a more symmetric form 



rX p 



(3.68) 



A 2 — [i 2 

As for the statement of existence; F exists if and only if exists, which is equivalent to the 
invertibility of the mentioned operator by the IIKS general theory. Q.E.D 

4 Tau functions and Fredholm determinants 

Slightly generalizing the definition in [2] (which is itself a generalization of the notion of isomon- 
odromic tau function introduced in the work of Jimbo-Miwa-Ueno [15, 13, 14]) we associate, to the 
space of deformations of the Ricmann-Hilbert problems 3.1 and 3.2 the two differentials below. 



Definition 4.1 We define the two forms over the space of deformations of Problem 3.2 and Problem 
3.1 



wJd) := I TrfSI^'JMiir 1 ) ^ (4.1) 

J 2l7T 
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and 

u F {d) ^\J^ (n^'JMM- 1 ) ^ (4.2) 

where d denotes any deformation of the jump matrices, ' the derivatives with respect to the spectral 
parameter and the integration is extended to all the contours where the jumps are supported, 7 + U7_. 



In the cases in which these two differential forms are closed it is defined, up to a constant, the 
corresponding tau function given by 91nr S / r = uj s / r (d). 

A particular case (of great interest) of deformations is when the jump matrices have the form 

M(A; s) = c T ^M Q (\)c~ T ^ , (4.3) 

and T(A) is a diagonal matrix depending on deformation parameters, while AIq(X) is assumed 
independent of them. A typical case is T(A) = J2k=o Tk^ k and the diagonal matrices T% are taken 
as deformation parameters. 

The relation between the Definition 4.1 and Fredholm determinants is elucidated in the following 
two theorems, stated side-by-side for comparison. 



Theorem 4.1 Given an operator K, as in Section 2 and the Riemann-Hilbert Problem 3.1 (with 
the same y{(J-)) we have the equality 



7+U7- 



31nr a = / Tr (EZ^LdM M' 1 ) ^- = <91ndet(Id 7+ - K 2 ) (4.4) 



Theorem 4.2 Given an operator JC as in Section 2 and the Riemann-Hilbert Problem 3.2 (with 
the same y{(J-)) we have the equality 

ainr r = i I Tr (TZ^dM M" 1 ) ^- = <91ndct(Id 7+ + JC) (4.5) 



Proof of Thm. 4.1. In [3] (see Theorem 2.1) it was proved (for the case of scalar operators, 
but the proof does not differ significantly as we see below) that 

d In t h = d In dct(Id 7+U7 _ - N) (4.6) 



20 



where the integral operator N, acting on L 2 ( r y + U 7_) is the one expressed in (3.43). On the other 
hand we have the identity 



det (Id 



7+U7- 



N) = det(Id 7+ - T o Q) = det(Id 7+ - K?) 



(4.7) 



where the first equality follows from 



det Id 



J-7 + U7_ 



G 
T 



= det Id 



det Id 



g 
t 


t Tog 



id 



7- 




det 
= det (ld 7+ 



Id 



7+ 



(4.8) 
(4.9) 



and the second is (3.47). The above computation is formal inasmuch as one would need to prove 
that all the operators involved are of trace-class. To see that we now prove that both J 7 , g are 
trace-class in L 2 (7+ U 7-,C p ). Recalling their definition (3.45) we augment the Hilbert space as 
M. := L 2 ( r y + U 7_,C P ) © L 2 (R, C r ), and extend trivially the definition of J 7 , g to the augmented 
space. This allows to represent them as the composition of two Hilbert-Schmidt operators (thus 
immediately implying the trace class property). Indeed (for example for IF) we have the identity 
below 



EZiQEkMx^ (0X 7+ (M) f d( ^(0x,_ (0 Mv)x 1+ (m) 



2m - 



(C-m) 



(4.10) 



which follows from Cauchy's residue theorem by closing the (~ integration either in the upper or in 
the lower half-plane. This realizes IF as the composition of two operators between the subspaces 
i 2 (7+,C p ) -> L 2 (R,C r ) -> i 2 (7_,C p ), each of which is Hilbert-Schmidt: 



m / idci- 



Trf^^)^!^) 



< 00 > 



7+ 



\M / |dc|- 



Trf£|(/i)^ 2 (/i) 



(4.11) 



For the sake of self-containedness we shall re-derive (4.6) below, let us denote by 3 = [A, B] 



the two block-columns of 2 (of sizes 2r x r) and by E 1 = 
r x 2r). Then 



C 
D 



the block rows of S 1 (of sizes 



J Tr( [A'_,B'_]dr«a+)dA-^ Tr 



7+U7- 
C 
D 



dA 



f Tr{KL l E'_dMM- x ) 

J 'V-i-U'V- 

[A'_,B'_]dW_ J dA = 



I Tr (DA'Sr) dA - f Tr (CB'Sf) dA 

J 7+ J 7 _ 



(4.12) 
(4.13) 
(4.14) 
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where we have used that C and A are analytic across 7+ and D,B analytic across 7_. On the 
other hand the jump relations imply 

(Wr 1 f A(/x)r(/x)d/x 
J 1+ fi — X 



A(A) 



l r 


r B( M )?( M )d M 


Orxr 


L,_ \x — x 



, B(A) 



rxr 
1, 



(4.15) 



and these identities can be differentiated on 7+ (for A) and 7_ (for B). We thus have 

- / Tr (DA'9r) dA — / Tr (CB'dr) dA = 

J 74- J 7- 



J dxfdfiTi (b{X) B ^^ dr(X^j dA + J dxj d^Tv (c(\) 



A(/x)r(/x) 
( M -A)2 



<9F(A) 



(4.16) 
(4.17) 



On the other hand these two terms exactly compute 



- Tr 7+ (11+- o clF) - Tr 7 _ (11-+ o 90) = -Tr 7+U7 _ ((ld 7+ u 7 _ + ^) o = 

= <91ndct(Id 7+U7 _ - N) 



(4.18) 



where 1Z is the resolvent of the operator with kernel N already used in (3.43) and 1Z has been 
decomposed as in (3.49) and F, Q defined in (3.45). Q.E.D 

Proof of Thm. 4.2 We start analyzing the l.h.s. of the equation and observing that 91nrr is 
equal to i<91nrH plus an additional term. Indeed, from T = R(X)E(X) with R defined in (3.25), we 
find 



r _ lr , = s _ ls , + E -1 R -1 R , E = S -1 H / 

so that 

2<9hiTr - <91nr H = / Tr (TZ^dM M' 1 ) 



2A 



S^(A)(l-aOS(A) 



(4.19) 



dA 



7+U7- 



2l7T 



Tr (EZ^LdMM- 1 ) 



7+U7- 



7+U7- 



1 

2A 



Tr SI A (A)aiS_(A) 



drx- Orxr 



dA 

2in 

dA 



J ^Tr(~: 1 (A)a 1 ~_(A)a+dr(A))dA + ^ ^ Tr (Sl 1 (A)a 1 3_(A)a_S?(A)) dA = 

= / ^-Tr(S: 1 (A)CT 1 S_(A)CT + ar(A))dA= f \ Tr (D(A)aj A(A)<9r(A)) dA (4.20) 
J 7 , A y A 



here we used the symmetries 



B(A) = <riA(— A), 



D(A) = C(-A)cti. 
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We want to identify this last integral; recall the notation 3 1 = ^ and that the jumps for H T 
imply for the column C T 



C T (A) 



7+ 



Moreover, by definition of inverses we have 

C(A)A(A) = l r = D(A)B(A). 

Consider 



^(A)A T (A) 



7+ 

D T (g)r T (Qdg N 
- + (£ + M)(A-0 



iff (A)A T (A)D r (fl£; 2 (0 Ef[(Z)dE 2 (vL) + dEf(QE 2 (fi) 
A-£ £ + M 

£i(A)A T (A)D T (0^(0 dEl(Z)E 2 (fi) 

7+ 



^(A)A T (A) 
/i + A 



_ ^(A)A T (A) {c T {x) _ C T { _^ dmi)+ j 

, ^ T (A)A T (A) T _ ^ pT ^ + /■ 
M + A J 7 



7+ 



A-e 
t t 



(4.21) 

(4.22) 

d£ = 
d£ = 



€ + M ^ - A 



D T (£)r T (0d£ I ^(/x) 



iff (A)A T (A)D T (0^ 2 (0 OEl (0E 2 (p) 



7+ 







£ff(A)A T (A)D T 




A-C 




iff(A)A T (A)D T 





d£ 
dC 
d? 



Ef {X)A T (X)aiT) T (ii)dE 2 (fi) Ef(X)dE 2 (^) , /" iff (A)A T (A)D T (£)£ 2 (0 diff^ M 



7+ 



M + A A + fi 

(note that we used (4.21)). Taking the trace we have to set fj, = X and integrate over 7+: the last 
term in (4.23) can then be simplified as well 

Tr (iff (A)A T (A)D T (£)£ 2 (O^iff (0^2 (A)) 



dC4.23) 



dA /d£ 

7+ ^7+ 



dA /d£ 

7+ ^7+ 



Tr 



dA /d£ r T (A)A T (A) 

7+ ^ 7+ 



1 / 1 



Tr 



(A -£)(£ + A) 
Tr (r T (X)A^(X)T> T (QE 2 (QdEl (£)) 
(A + A) 

D T (0£ 2 (£)a^(e) 
D -r (03»(0^a(0d£ 



(4.24) 



2^VA-C A + C 
B r (0-B r (-0„ 



7+ 



(4.25) 
(4.26) 
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Tr ((l r - A T (0SiT> T (H))E 2 ($dE?($) 



7+ 



2i 



(4.27) 



Taking the trace of (4.23) we thus have 



=Tr dlC 



Tr (ft ++ o dK) = - J dA * — '- 



Tr 



A T (0^iD r (0^(0^(0 



7+ 



2£ 



7+ 



A T (Qg 1 D T (Qi; 2 (c)^r(c) 

2£ 



Together we thus have 



J 7+ 2A 



and therefore 



Tt(YZ x V_dMM- 



7+U7- 



' 2m 



Tr (2I 1 2'_9MM 



7+U7- 



dA 

J 2i7T 



(4.28) 
(4.29) 

(4.30) 

(4.31) 
(4.32) 



= 2 Tr(7e++ o dK) + 2 Tr dJC 
In summary, using Theorem 4.1, we have 

/ Tr (TZ^dMM- 1 ) — = dlndct(Id 7+ - K 2 ) + 2Tr(ft ++ o dK) + 2TrdK (4.33) 



On the other hand we now show that the r.h.s. of (4.33) is precisely 291ndct(Id 7+ + K) at which 
point the proof shall be then complete. To verify this last point we have (using Theorem 3.1) 



(Idy. +K)- 1 



= (Id 7+ - /C 2 )- x (Id 7+ -/C) Th = 3I (W 7+ +^ ++ )(Id 7+ - K) 
= Id 7+ — K + — 1Z++K 



(4.34) 



from which we can compute the variations of the determinant 

2<91ndet(Id 7+ + K)= 2Tr((Id 7+ +K.y 1 dK) = 

= 2 Tr {dK - KdK + TZ++8K - TZ++K8K) = 

= - Tr((Id 7+ + TZ++)d{K 2 )) + 2 Tr(dK) + 2 Tr(TZ++dK) = 

= <91ndct(Id 7+ - K 2 ) +2Tr(n ++ dK) + 2Tr(<9/C) 



(4.35) 



Q.E.D 

When the dependence on the deformation parameter is in the form specified in (4.3) then we 
can write the differentials in terms of formal residues as shown here. 
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Proposition 4.1 Suppose that M(X) in (3.2) can be written as M(A) = e T ' A ^Mo(A)e _T ^ A ^ where 
T(A) is a polynomial diagonal matrix without constant term in A (T(0) = 0) whose entries depend on 
the deformations and such that T(— A) = a\T{X)a\. Let d be the derivative w.r.t. one deformation 
parameter. Then 



where the residues are understood as formal residues, or the coefficient of X 1 in the expansion at 
infinity. 

Proof of Prop. 4.1. The equivalence of the formal residues (4.36) with the integral represen- 
tation (4.1) (or (4.2)) was proven in [2] in a more general context, but we recall here the gist of it. 
The formal residue in (4.36) (for the case of Wr, the other case being completely analogous) can be 
written as an integral on an expanding counterclockwise circle (with the piccewise-defined T) and 
then it can be transferred by the use of Cauchy theorem to the integral 



Firstly, the term crossed out is zero because T(A) = 0(A) => <9T(A) = 0(A) (as A -> 0) and T^V 
may have at most a simple pole at A = (thanks to (3.10)) so that the term is analytic at A = 0. 
Secondly, note that M^ 1 ^!' is (piecewise) strictly upper or lower triangular on 7+ U 7_ and dT 
is diagonal, hence the term Tr(M -1 M'<9T) = on 7+ U 7_. On the other hand from the formula 
(4.3) follows immediately that MdTAI^ 1 - dT = —dMM~ x and hence (4.36) gives exactly (4.2). 
Q.E.D 

The following corollary follows from direct matrix multiplications using the asymptotic forms 
(3.9) and (3.7) together with the special structure of the expansion matrices (3.15) and (3.21). 

Corollary 4.1 If dT(X) = iXekk®o~3 with the diagonal elementary (r x r) matrix, then 




(4.36) 



(4.37) 




(4.38) 



(4.39) 



(4.40) 




(4.41) 



25 



Similarly 

0Js{d) = - res ( E^E' \e k k®o- 3 )dX = -2i(oti)k,k (4.42) 



To conclude this section we prove that, when r = 1 the relation between the two Riemann- 
Hilbert problems, and in particular equation (3.23), can be interpreted as a Miura transformation 
between the two tau functions. 

Proposition 4.2 Suppose r = 1 (but possibly p > 1) and M (A) := e isX " 3 M (X)c- isXa3 (i.e. a 
special case of Prop. 4-1). Then the tau junctions (Fredholm determinants) r=,Tr are related 
through the Miura transformation 

(d s In t s - 2d s In t v ) 2 = ~d 2 s In r s (4.43) 
Equivalently we may simply write 

(wa(0.) - ^ r (d s )) 2 = -0 s ^(d s ) (4.44) 

Proof of Prop. 4.2. The solution 2 of Problem 3.1 is such that 

d s (S(A; s)e lsA<T3 ) =U{\; s)E{\; s)c lsXs3 (4.45) 
U(X; s) := iAcr 3 + 2^ (s)oi (4.46) 

which can be easily proved by noticing that C/ has no jumps, hence it is entire, and then by looking 
at the behavior at infinity using the expansion of E. On the other hand then comparing the terms 
in the expansion of the two sides of (4.45) one finds that 

d sCn = -2i(3 2 . (4.47) 

Now, Corollary 4.1 (i.e. Prop. 4.1) yields 

a s lnr = = - res Tr(E- 1 E'd s T)d\= -2ia! (4.48) 

" A— oo 

<9 s hiT r = -- res Tr (r _1 r'9 s T) dA = -ia x (4.49) 

2 A— oo 

Rewriting (3.23) as 2fi\ = {2ia x - 2iai) 2 and using the equations (4.47), (4.48), (4.49) we obtain 
the statement of the proposition. Q.E.D 

Remark 4.1 Defining u := 2d 2 In if and v 2 := —d 2 liiTs we obtain the usual formulation of the 
Miura transformation 

u = —v 2 ± d s v. 
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5 Applications: Fredholm determinants and noncommuta- 
tive Painleve II, XXXIV 



Wc now consider the Fredholm determinant for the convolution operator on L 2 (M + , C r ) given by 



Ai(x + y;s)f(y)dy 



Ai(x;s):= [ e ^Ce e ^e ix ^ = [c jk M(x ■ 
J 1+ 2tt 



s k)} 



1, 



ISljJL 



lS 2 fl 



ISrfJ- 



ISfl 



s := diag(si,s 2 , ...,s r ) 



(5.1) 
(5.2) 

(5.3) 
(5.4) 



The matrix C is a constant r x r matrix and the contour 7+ is a contour contained in the upper 
half plane and extending to infinity along the directions arg(z) = ^, ■^ L . 
Here the matrices E\ , E2 can be chosen as 



£i(A) = -^-e e ^C, E 2 (X) = e 9 ^ , r(A) = 



1 



(5.5) 



The first issue is whether the solutions of Problems 3.2, 3.1 exist for real values of the parameters 
s. We shall -in fact- show an existence theorem for Problem 3.1, which immediately implies existence 
of the solution of Problem 3.2 by Proposition 3.2. 



Theorem 5.1 Suppose C — is a Hermitean matrix; then the solution to Problem 3.1 with r as 
in (5.5) exists for all values of s G R r if and only if the eigenvalues of C are all in the interval 
[—1, 1]. If C is an arbitrary complex matrix with singular values in [0, 1] then the solution still exists 
for all s€l r . (The singular values of a matrix are the square roots of the eigenvalues of C^C) 



Proof of Thm. 5.1 The proof is based on the estimate of the operatorial norm of the operator 
Aig (with parameter C G Mat(r x r, C)); the invertibility of the operator Id+^li^ will be guaranteed 
if the norm of Aig is less than one. 

On the other hand the invertibility is equivalent to the non- vanishing of the respective Fredholm 
determinant; hence from Corollary 2.1 we have (in the present notation, with p = r) 

det (ld R+ ± Ai s ) i2(a+>cr) = det (Id 7+ ± K) LHl+ ^ P) (5.6) 
det (ld R+ - Aii) L2(K+iCr) = det (Id 7+ - K 2 ) i2(7+jCP) (5.7) 
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Thus, if || 1.4^11 1 < 1 then |||-4if||| < 1 and thus the Fredholm determinants on the line (5.7) do not 
vanish. This is sufficient for the existence of the solution of Problem 3.1 as shown in Thm. 3.1. 

Let us then estimate the norm |||4is|||; first of all note that L 2 (K+) ~ L 2 ([s, oo)) by simple 
translation; with this in mind we can express the operator Ais as the operator Aiq but acting on 
the space 

Hg= L 2 {[ Sll oo)) © . . . 8 L 2 {[s r , oo)) (5.8) 
Ai 5 -.Us -> Us (5.9) 

(/i,---/r)^ \ f2c jk ^Ai(x + y) Xl s h ,o )fk(y)dy\ (5.10) 



\k=l 



Let Pg be the orthogonal projector 

P?:L 2 (RX r )^Hs, P ? =diag( X[sii00) ,...,X [Sr , oo) ) (5.11) 

Then we have Ais — PsAi^Ps- On the other hand it is evident that the operator Aiq : L 2 (M., C) — >: 
L 2 (M, C) is the tensor product Ai$ = C<E)C where we have denote C the scalar convolution operator 
with the Airy function on R + 

C :L 2 (E) -> L 2 (K) (5.12) 
£/(x) := / Aii.r • / /:/: / /;d // (5.13) 



This operator squares to the identity (as it is easily seen in Fourier transform, but is also well known 
[22]) and hence has unit norm (in fact it is a unitary operator, an easily verified fact in Fourier 
transform) . Therefore 

|||4%||| = |||C®£||| = |||C||||||£||| = |||C||| , (5.14) 

and hence 

||L4z 3 ||| = |||PM%Pr||| < Ml HIAjIII W\Ps\\\ < IIICHI (5.15) 

Since |||C||| is the maximal singular value the first part of the theorem is proved because |||C||| < 1 
implies that the norm of our operator is strictly less than one. 

To prove necessity in the case C = Hermitean, suppose that C has an eigenvalue >c € 
K \ [—1,1] with eigenvector vo (for Hermitean (and more generally normal) matrices the singular 
values are simply the absolute values of the eigenvalues). We need to show that for some choice of s 
the Fredholm determinant vanishes; we will accomplish this by finding a special value of s for which 



28 



the square of Aig has an eigenvector and hence it is not invertible, thus implying the non-solubility 
of Prob. 3.1. To this end we take s = (s, s, s, s . . . , s) and f(y) = i> Q ip(y) with ip(y) <G L 2 (R + ). Then 

(Aijf)(x) = k 2 v / K Ai {x + s,y + s)(p(y)dy (5.16) 

where K A i is the well known Airy kernel 

T , . , Ai(x)Ai'(y) - Ai'(x)Ai(w) /" ... , ... , , 

K Ai (x, y) = ——^ — — — — = / Aifx + z)M(z + y)dz (5.17) 

x-y J R+ 

It is well-known that K Al on L 2 ([s, oo)) is self-adjoint and of trace-class. Let A(s) be the maximum 
eigenvalue. This is a continuous function of s and tends to 1 as s — > — oo (it clearly tends to zero as 
s — > +oo) [22]. Let f s (y) be the corresponding eigenfunction and use now f(y) = VQ(p s (y). Then 

(Aiif)(x)=A(s)K 2 f(x) (5.18) 

If k 2 > 1 there is a value of So S R for which (^li|/)(x) = f(x), thus proving that Id — Ai 2 cannot 
be invertible for s= (sq, . . . , sq). This concludes the proof. Q.E.D 

The reader may now wonder whether these kernels have a "physical" interpretation. The answer 
is in the affirmative 

Theorem 5.2 IfCis real or Hermitean then the kernel Ai 2 ? is totally positive on the set {1,2, . . . r}x 



Proof It is a general result that if (X, d/i) is a measure space then for any W C X we have 

fcldct / / a (C)fl6(C)dMC) = / dct[f a (C c )]dct[g b (C d )]f[dfi(C 3 ) (5.19) 

In our case we take A = {l,...,r} x R with the counting measure times Lebesgue measure. 
A function on X is then equivalently interpreted as a vector of usual functions: £ = (j, x) =>■ 
/(C) = /((i) x )) =: fj( x )- With this understanding the kernel Ai| is understood as a scalar function 
on A x A to wit 

•2 



^(6,6)= (5.20) 



[Aig]j 1 ,j 2 (xi,X2) =y2c jlk c k j 2 / d^Ai(a;i + z + + s k )Ai(x 2 + z + s 3 - 2 + s k ) = (5.21) 
fc=i Jr + 

= f MQFfoiQF&b) (5-22) 
Jx + 
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where we have set 



'■= c jltk Ai(xi + s n + z + s k ) 



( = (k,z)eX+ :={!,. ..,r} 



C X 



The statement of total positivity amounts to checking that for any K <E N and any £1 , 
we have 

det |y«Kf«,&)] 0)i < JC >0 
To this end we use the previous fact (5.19) and we have 



det [Aig(£ a , &)] 1<a 6<K = det 



d^(e)F(^,C)^(C,6) 



(5.23) 
(5.24) 

• •,6r ex 

(5.25) 
= (5.26) 



K 



K 



= Y\I dct[F(^,Cc)]det[F(C c ,Ca)]n d M(Cc) = ^y / Jdet [F&, ( c )}\ 2 f[ dfi(Q > (5.27) 

where the modulus occurs if C is complex Hermitean (in which case F(£,£) = ^(C>£))' while if C 
is any real matrix then we have a simple square (which is anyway positive). Q.E.D 

Theorem 5.2 allows us to interpret the kernel Aig as defining a determinantal point process 

on the space of configurations X = {l,...,r} x K [21]. The Fredholm determinant is then a 
multi-level gap distribution for said process on the interval [S, oo) (after a translation x i-> x — S). 

5.1 Noncommutative Painleve II and its pole-free solutions 

We consider first Problem 3.1 for 3; the jump is written (r defined in 5.5) 

5+ = 3_ (l 2r - 2i7rr(A)»cr+) , A 6 7+ 
3 + = 3_ (l2r — 2nrr(— A)®<7_) , A e 7- 
The matrix *(A) := 5(A)e 0(A)<SCT3 , with 6»(A) as in (5.3), solves a RHP with constant jumps 

f + = (l 2 r + C«CT+) , A € 7+ 

*+= (l 2r + C®ct_) , AG7_ 
*(A)= (l 2r + 0(A- 1 ))e e(A)8<T3 , A^oo (5.30) 

It would be simple to show that ^(A) solves a polynomial ODE in A (of degree 2, see Lemma 5.1), 
which eventually would lead to showing that /?i(s) solves a noncommutative version of the Painleve 
II equation (whence the title of the section). In this perspective, the above jumps are a particular 
choice of Stokes' multipliers associated to such an ODE, exactly as in the scalar commutative Lax 
representation of PII [12]. We thus describe in the next section below, ex ante, the most general 
set of generalized monodromy data for the ODE (5.46). 



(5.28) 
(5.29) 
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5.1.1 The general Stokes' data/Riemann Hilbert problem for 

Denote by A C C r the set of diagonals 

A:={seC r : 8j=8 t , j^i} ■ (5.31) 

Let s^°) G C r \ A and choose a ray 7^ := M + e~ Vi? in such a way that ^(sj. ' — s^ ') 7^ for z G 7^. 
Let 7l := — 7h- We introduce the ordering k -< ^ as follows 

k<£ ^ S(e WB (4 0) -4 0) )) <° ( 5 - 32 ) 

For a fixed 7^ this ordering is constant in a suitable open conical neighborhood of s*- -* not 
intersecting the diagonals A (as should be clear by a simple continuity argument): we shall under- 
stand such choice of neighborhood and keep the chosen ordering fixed. We shall say that a matrix 
TV is upper(lower)-triangular relative to the ordering (5.32) if N^e = for k -< £ (£ -< k, 
respectively) and = 0. 

Example 5.1 If s\ < S2 < ■ ■ ■ s r are real and ordered, then the notion of upper (lower) triangularity 
relative to any ray argz = (— tt, 0) is the usual one. 

We define the six additional contours 

:= I + eT+f , fc = 0,...,5 (5.33) 

Let Co, . . -C-2 three arbitrary r x r matrices, S u = l r + N u , Si = l r + TV; with N u ,Ni two 
upper/lower triangular matrices relative to the ordering (5.32) determined by the choice of jr, and 
let M = diag(/ii, . . . , fi r ) G SL(r, C) (traceless). The entries of M will be referred to as exponents 
of formal monodromy. 

Problem 5.1 Let ^f(X) be a sectionally analytic function on C \ 70 U . . . 75 U 7l U 7r, bounded over 
compact sets of C and solving the following Riemann- Hilbert problem 



*+= 




-(l2r- 


I- C ®CT+) , 


A G 70 


*+ = 




-(l2r 


+ Cl®(7_ 


-) , 


A G 


Ti 


(5.34) 


*+= 




-(l2r- 


f- C 2 S)Cr + ) , 


A G 72 


*+ = 




-(lar 


+ Co®CT_ 


-) . 


A G 


73 


(5.35) 


*+= 




-(l2r- 


)-Ci8CT + ) , 


A G 74 


*+ = 




-(lar 


+ C2®(T_ 


-) , 


A G 


75 


(5.36) 


*+= 




-[(lr- 




-S)] ■ 


A G 7_r 














(5.37) 


*+= 




-[(lr- 




£.)] , 


A G 7l 














(5.38) 


*+= 




_ e - 4 ^M®l 2 ^ A g 


R± 
















(5.39) 


*(A)= 


(1 


2r + C(A -1 )) e -^M8l JA M8X Je )W 8aS 




A -> 


X) 








(5.40) 
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where in (5.40) e = 1 in the upper half-plane, e = in the lower half-plane and arg(A) G [— 7r, it). 
The matrices Co, . . . , C r , S u , Si, M are chosen satisfy the no-monodromy condition stating that 
the product of the jumps is the identity. (We choose 7# 7^ R± and all the rays are oriented towards 
infinity). 

Since the rays , 1r may lie in between different jj 's depending on the value of s the no-monodromy 
condition may take different forms. For example, ifs£l r \Awe can choose arg(7#) = ^ + e and 
the no-monodromy condition takes the form 

(l ar + Co®CT + )(l 2r + C 2 «o-_)(l 2r + C 1 c Sf r + )e- raM x 

(l ar + C ®<7_)(l 2r + CW+)(l 2r + Ci« . <r_ )e~ ilrM = l 2r (5.41) 

Remark 5.1 XTie problem associated to the Fredholm determinant of the operator as in Thm. 5.1 
corresponds to the particular choice S u = Si = C\ = M = and Cq = C = — C 2 . 

Note that the jumps satisfy the symmetry M(— A) = diM(X)&i and hence we also have (noticing 
that 9(— A)®i73 = ai6(\)<sa^5i since 0(— A) = —0(A) as per (5.3)) 

*(-A) = cti*(A)cti (5.42) 

The dimension of the manifold (Co, C\, C 2 , S u , Si, M) of solutions of (5.41) can be computed by 
noticing that there are a total of 3?* 2 + 2 r ^~ 1 " 1 + r — 1 = 47- 2 — 1 variables. The equation (5.41) is 
of the form An\A5\ = l 2r and hence -due to the symmetry of conjugation by a±- there are only 
2r 2 independent equations. Of these, one is redundant since the determinant of A is already unit. 
Hence there are 2r 2 — 1 independent equations and thus the manifold of solutions has dimension 
2r 2 . 

Lemma 5.1 Let the matrix M/(A) be the solution of the Problem 5.1 and denote the asymptotic 
expansion at 00 as 

00 „ 

^ A ) e ie*M 9 l aA -M»l Se -*(A)8>os = ^ + 2l , 

H2J+1 = a2 3 +i®(T3 + /3 2 j+ic 2 , S 2 j = a 2j ®l + /3 2 jCi , (5.43) 

(recall that e = 1 /or QA > and e = /or QA < 0) where the expansion is valid sectorially and 
independent of the sector. Then 

d Sj ^= Uj^ , Uj = i\ej®s 3 +i[ai,ej]®l + {0i, ej}®<ri (5.44) 

A r 1 

d\$>= A(X)H> , A(X) = - ^ Uj - -D (a v »a 3 + f3 v »a 2 ) + is®cr 3 = (5.45) 

" 3=1 
A 2 1 

A(X) = i—a 3 + Aftsai - -D/3i®<7 2 + + s)®a 3 (5.46) 



l + 5„ 
1 + Si 



1 + Si 

1 + 5., 
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where 

r 

B :=J2 d S] , ej := diag(0, 0, . . . , 1, 0, . . .) (5.47) 

i=i 

with the one in the j-th position. 

Proof The fact that the expansion for (recall that arg(A) £ [— tt, 7r)) 

f ^(A)e i7rM A- M55l2 c- e ( A )®^ 3A > , |A| > 1 
5(A) := < ^(A)A- Ms l2 c- e ( A )® CT 3 QA < , |A| > 1 (5.48) 
I *(A) |A| < 1 

near A = oo is of the form in (5.43) follows from the symmetry \&(— A) = ai^(X)ai which then 
implies the same symmetry for 5. The function 5 has then no jumps on K± \ {|A| < 1} and the 
remaining jumps are those of ^ conjugated by e -leM A M ® l2 e 6 '( A )®' :r3 . The fact that the expansion 
is independent of the sector is a consequence of the fact that the jumps for S along the eight rays 
are analytic in a small open sector around said rays and of the form 

S + (A) = S_(A)(l + 0(A-°°)) , A^oc, Ae 7 oU...75U7 L U 7i? (5.49) 

uniformly within said sectors. The fact that Uj and A are polynomials is an immediate consequence 
of the fact that the jumps of ^ are independent of A, s. Using Liouvillc's theorem and the fact 
that 9 Si >I n I'~ 1 is entire (a simple consequence of the independence on sj of the jumps) we deduce 
immediately that Uj(z) can only be a polynomial of degree 1. Then 



3 -S(A) + /AE; Air, a 3 = {uf ] \ + t/j 0) ) 5(A) 



<9,(ai®<7 3 + /3i®<72) ., /„ ai»<73 + Bi®a\ «2®1 + /32»ci 

+ l \ 1 + + + . . . ej ®a 3 



A V A A2 

= ([/j 1} A + C/f) (l + '" "•■ A " + " 28l +/ 28fTl + . . .) (5.50) 

Comparing the coefficients of the powers of A we have 

A : => Uf> = i\ej»a 3 (5.51) 

A : => i [ai®er 3 + /3i®02, Cj»er 3 ] — £/j°' (5.52) 

A -1 : => d Sj (ai»a 3 + /3i»cr 3 ) = —i [a 2 »l + ej®cr 3 ] + (ai«>er 3 + /?i®cti) 

d S] (ai9<7 3 + f3i<sa 3 ) = —i [aa®l + /?2®cr2, ej®cr 3 ] + uf^ (ai®cr 3 + Pvsa 2 ) (5.53) 

If we sum up for j = 1, . . . , r we obtain the differential equation 

D (ays,a 3 + /3i®CT2)= — i [«2®1 + /?2®cr2, 1®o- 3 ] + i [ax®os, +/3i®er 2 , l®cr 3 ] (ai®cr 3 + /3i®cr 2 ) = 

= -2/3 2 »cri+2^iscri (ai®cr 3 + ^i»CT2) = — 2/3 2 ®cri + 2i/3iai®c72 — 2if3f»a 3 (5.54) 
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In particular 

Dai = (5.55) 
If we look also at the A~ 2 coefficient we find 

Dft = 2z/3 2 - 2(3 iai , (3 2 = - ^D/3i - iftai . (5.56) 
Exactly as before we argue that A(z) is a polynomial of degree 2. Then we compute 

9 A S(A) + 3(A) (iy 1 P + is) 8(73 - {A 2 X 2 + A x A + A ) 5(A) (5.57) 

/ 0-1003 + /3i®cr 2 a 2 '»l + A2®0-:l . \ /.A 2 . \ , . 

...+ 11 + + ^ + ... I Uyl r +«Sj *(T 3 = (5.58) 

= [A 2 z +Aiz + A )\l-\ 1 h...J (5.59) 

Collecting the coefficients 

A 2 : =*> A 2 = ^X 2 Ua 3 (5.60) 

A 1 : =>■ A\ = -[o!i8£73 + Pi»(J 2 , Ibcts] (5.61) 

A : =>■ Ao = is t :) + -[a 2 ®l + /?2®0i, l«cr 3 ] — Ai (a<ig><73 + /3i®0"2) = (5.62) 

A = is«(7 3 - -D (ai®a 3 + f3i»cr 2 ) (5.63) 

where we have used formula (5.54). The second expression for A(X) (5.46) follows from (5.55). The 
rest of the proof is a simple computation. Q.E.D 

Lemma 5.2 Let ^ be as in Lemma 5.1 and denote by j3\ = Pi(s) the r x r coefficient matrix 
in Si = ai0 03 + /3i®<72 of the expansion as in the mentioned Lemma. Then the matrix function 
/?i(s) £ Mat(r x r, C) satisfies the noncommutative Painleve LL equation. 

r q 

D 2 /3i =4{s,/3i} + 8/3? , s:=diag(«i,...,« r ) , D := ^ - — , = + (5.64) 

3=1 ^ 

Proof. We use the zero curvature equations 

r 

D* = ^ t/j* =: 17 D * , <9 A * = A* (5.65) 

3=1 

(5 A C7 D + [/ D A-DA + AC/ D ) = (5.66) 



34 



with the Uj's introduced in Lemma 5.1. We have 

Ur> = 'Al er 3 + 2/3x®<ti , A = -jU-& — — D {a\mu^ + /3i®er 2 ) + js<»cr 3 (5.67) 

j 1 

A = -A 2 l»cr 3 + \f3\®ai - -D (ai®cr 3 + /3i<8cr 2 ) + is»a 3 (5.68) 

We now compute this expression (for simplicity we denote with a prime the action of D. noting 
that Ds = l r ) 

d\Ur>= «l®cr 3 (5.69) 

DA= - i (a'/^crs + /?" »<j 2 ) + «W 3 (5.70) 

i\l®a 3 + 2(3 1 «cr 1 , — - (ai®(T3 + /?i©CT 2 ) + zs®cr 3 



(5.71) 



A/3'^tri - 2ij3iai®cr2 +t{^i,j8i}® cr 3 - 2{^ ls s}® cr 2 (5.72) 



Hence 



0=<9 A £/ D -DA+ [f/ D ,A] = (5.73) 
=-a' 1 '®a 3 + -jflJWa - 2i / 1 a / 1 ®cr 2 + /3i}®<7 3 - 2{/3 l5 s}®o- 2 (5.74) 

Using now (5.55) we have a" = —2i{/3i, /3[} and hence we are left only with 

- 4/3? - 2{ft, s}) 8 a 2 = (5.75) 

Q.E.D 

Thus Lemma 5.1 and the matrices (5.65) provide a Lax matrix representation for the general 
solution of the noncommutativc Painleve equation (5.88) which is parametrized by the 2r 2 initial 
values Pi(s) and D/3i(s) at any point s. 

It should also be clear that any solution j3\ of the noncommutative Painleve II equation (5.88) 
yields a compatible Lax pair A(X; s) (5.46) and ?7d(A; s) (5.65); the Stokes' phenomenon (generalized 
monodromy data) for the ODE d\fy = A^> can be seen to be given exactly by the data specified 
in Problem 5.1 (we refer to [24] for the general theory of Stokes' multipliers). Thus any solution of 
(5.88) is obtained via the above Lax-pair. 

Remark 5.2 The generic solution of Pi(s) of the noncommutative PII equation will have non- 
movable singularities on the diagonals s € A; this is due to the presence -in general- of nontrivial 
Stokes multipliers along the rays 7i?,7L- If those multiplier are trivial as well as the exponents of 
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formal monodromy i.e. M = 0, S u = = Si then those singularities will be absent. In this case the 
no-monodromy condition (5.^1) can be spelt out more clearly as 



[C ,C 1 ] = 0, [C 1 ,C 2 ] = 0, [C ,C 2 ] = 



(5.76) 
(5.77) 



C + C 2 + d + CodC'2 = 



which yields a manifold of dimension r 2 + r ( the matrices can be generically diagonalized simulta- 
neously and a simple counting yields this number). The condition (5.77) resembles very closely the 
ordinary situation r = 1 of the commutative Lax representation for PII [10] 12 . 

The importance of the isomonodromic representation for the noncommutative Painleve II equa- 
tion is that it implies automatically the Painleve property [18] that the only singularities of the 
solution are poles except -possibly- the singularities on the diagonal manifold s 6 A if the Stokes' 
matrices S u ,Si are nonzero. 

Remark 5.3 Another important remark is that the solution /3i(s) as a function of the barycentric 
variables 



has only poles as a function of S (note that D = &s) if s ^ A; this is so because changing S does 
not change the differences between the Sj 's and hence never crosses the diagonal manifold. 

Remark 5.4 To our knowledge, the noncommutative Painleve equation (5.64) has appeared first 
in the recent [19] where the authors construct special rational solution using the theory of quasi- 
determinants [7]. Previously, a version with scalar independent variable (hence replacing the anti- 
commutator by simply sf5\) was studied in [1], where the Painleve test was applied. It seems that 
the Lax representation for the noncommutative version of [19] appears in the present manuscript 
for the first time. It seems possible to generalize the Lax-pair representation by allowing a pole at 
X = in the Lax matrix A(X) (exactly as in the scalar case). For example the compatibility of the 
following two Lax matrices 



Their matrix ^(A) has the symmetry A) = ct2 , J'(A)(T2 , which means that it should be compared with ours 




(5.78) 



A(X)= —d 3 + A/?i®eri + i (s + (3 2 ) »cr 3 - -P[»a 2 + -6cti 
[d x -A(X),T)-Un(X)} =0 



Ur> = iXa 3 + 2[3 V : i ,a 1 (5.79) 
(5.80) 



after conjugation by e 1 4 
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with & an arbitrary scalar (i.e. commutative symbol). The zero- curvature equations are easily 
verified to yield 



D 2 /?i = 4{s, Pi} + 8j3f — 46 



(5.81) 



which is precisely (with different symbols) the Painleve II equation studied in [19]. From the isomon- 
odromic method, however, the above equation appears to be not the most general that one may obtain 
by allowing a pole in A{\). We do now dwell further into the matter since it is peripheral to the 
focus of the present paper. 

The compatibility equations for the operators d Sj — Uj, d Sk — U k and d\ — A yield additional 
equation listed in the Corollary below, which is proved along the same lines (but we will not report 
the proof here since it is unnecessarily long, straightforward and anyway this has no bearing for our 
goals) 

Corollary 5.1 The matrices /3i,a\ satisfy the systems (dj = d Sj , D = ^2 dj, j = 1, . . . , r) 

1 



djP\ =-^{^o, D/?i} - i{ej,piai} + iej[ai,fii] + ioaejfii + iPiejCti 
1 i 

-djDPi =i(dj (ii)ai + {{/3 1} e-,}, s} - -[ej,ax]Dj3 x + [[ati,ej],Pi]ai ^ 



1 - Cj + i 



s, [ai,ej 



(5.82) 

(5.83) 
(5.84) 



{ej,dkPi} - {ekjdjPt} =i[p x e k , ai&j] +i[ej/3 1 ,e k ai] + ie k [0i,ai]ej + 

+ie J [ai,/3 1 }e k + i[aie k ,Piej] + i[ejai,e k Pi] , (5.85) 
i[e kl dja{\ - i[ej,d k ai] ^kafcj - Cja\c k + e k /3fcj - Cj/3fc k + 

+ [e k /3 1 ,e j l3 1 } + [/?ie fc , frej] + [e^-a-i, e k a{\ + [a^j, aie k ] . (5.86) 

5.2 Pole-free solutions of noncommutative Painleve II and Fredholm de- 
terminants 

We now return to the specific situation of the RHP associated to the integrablc kernel Aii; this is 
the special case of the setting as explained in Remark 5.1. 
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Theorem 5.3 Let S = S(A; s) be the solution of Problem 3.1 with r as in (5.5); let 

/3i(S) := -» lim AH 12 (A;s) (5.87) 

A— >oo 

where denote the r x r blocks of 3, i,j = 1,2. 27ie matrix function /3i(s) € Mat(r x r. C) 
satisfies the noncommutative Painleve II equation. 

r n 

D 2 /3x =4s/3 1 +4/3 1 s + 8/3 1 3 , s := diag(*i, . . . , « r ) , D := ^ — (5.88) 

j=l 3 



The asymptotic behavior of the particular solution associated to the Problem 3.1 is as follows: if 

1 yvr 



S := - s j ~^ +°° an d &j '■— s j ~ S , j = 1, . . . ,r are kept fixed, \8j \ < m, then 



[fain = -c ke Ai(s k + s£) + LfSe-^ 2S - 2m ^\ (5.89) 

If C is Hermitean then so is the solution P\(s) of the noncommutative Painleve equation (5.88) and 
it is pole- free for all s£l r if and only if the eigenvalues of C are within [—1,1]. If C is arbitrary 
and its singular values lie in [0, 1] then the solution is also pole free for s € M. r . 
Finally, the Fredholm determinant r(s) := det (id — Aig 2 ) satisfies 



dct (Id - Ai/) = exp ^-4 J (t-S) Tr(f3 2 (t + <f))dij 
where t + d := (t + Si, ■ ■ ■ , t + S r ). 



(5.90) 



The last statement is the noncommutative (matrix) equivalent of the celebrated Tracy- Widom 
distribution [22]. Before giving the proof of Thm. 5.3 we prove the uniqueness of the solution. 

Proposition 5.1 For any r x r matrix C there is a unique solution of noncommutative PII (5.88) 
with the asymptotics (5.89). 

Proof of Prop. 5.1. The proof docs not differ significantly from the scalar case as in [9]. In 
baryccntric and relative coordinates S, 5j as in Thm. 5.3 we have D = d,s- The regime we consider 
is S — > +oo and all Sj bounded below. We note that the function 

[U(s)} ke := -c M Ai(s k + s e ) - -c ke Ai(2S + S k + 6 e ) (5.91) 

is a solution of the linear part of (5.88): 

T> 2 U k e = -4c ke Ai"(s k + s e ) = -4c k ,t(s k + s e )Ai{s k + s e ) = A(s k U kt + U ke s e ). (5.92) 
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Then any solution {j3\)ki with the specified asymptotic also solves the integral equation 13 

poo 

[Pi)kt = U ke +4n (Ai(2S+S k +Si)Bi(2t+5 k +5e)-M(2t+6 k +5t)Bi(2S+5k+5i)) pi) 3 ] fef d<5.93) 
Js 

Equation (5.93) can be solved by iterations for S sufficiently large, as noted in [9] for the scalar 
case. 

The local uniqueness follows from the local uniqueness of the solution of the ODE (5.88) (in S). 
The solution is easily seen to be locally analytic in s because of the analyticity of the ODE and also 
from the integral equation. We also point out that since the generalized monodromy data associated 
to this solution have S u = Si = = M (see Problem 5.1) then there are no critical singularities at 
all in s and we may only have poles at most (the Frcdholm determinant is analytic in s, hence can 
only have zeroes). Thus the solution is globally defined for s E R by analytic continuation. Q.E.D 

Remark 5.5 Because of (5.89) and Prop. 5.1 we may call the special solution of noncommutative 
P II arising above the noncommutative Hastings-McLeod solution(s). 

Proof of Thm. 5.3. The fact that f3i solves the noncommutative PII equation (5.88) follows from 
Lemma 5.2 since this is a special case of that with S u = Si = C\ = M = and C = Cq = —Ci- 
Asymptotics. Suppose that S — £ ^2 s j is large and positive and 5j := Sj — S are bounded by 
-let's say- m. We rewrite the RHP in the scaled variable z := -^=. The jump on the contours j± 
of the form 1 — 2kr« er+ , 1 — 27rf » er_ can be factored into (commuting) matrices (here below e k t is 
the elementary matrix) 

1 2/-r ,t. f[ (l + c M e S Hi z3+{2+ ^ 1)z )e k jc i <T + ) (5.94) 

kj=l ^ ' 

1 2/xr ,7 f[ fl + c H e" ,si (' z3+(2+ ^ )2 )e fc ^(T-) (5.95) 



k,i=l 



Each factor has a saddle point at z = ±iy2+ Sk + Se and the contours j± supporting the single 
jump can be split according to the factorization (5.94) so that each of the factor is supported on a 

(k £) 

different contour 7± of steepest descent for the corresponding phases. Proceeding this way the 
reader realizes that each factor 

<l ' (v Sz) = 1 + ci-pc V 3 V s J J ei. e(xo-+ . z e Vl ' 



M^'(VSz) = l + c ke c^^y* z ^—JV ek , m a ± , z€ 7 f (5.96) 
is close to the identity jump in any L p (j±'^), 1 < p < oo with the supremum norm given by 

||1 - A/ (M) ||oo = | C M \ C -I^+S k+ S e )i _^ _ (5 9?) 



13 We recall that Bi is the solution of /" = xf such that Wr(Ai, Bi) = it 1 . 
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This shows that the RH problem may be solved by iterations; the first iteration for S yields 
J 1+ (w-z) } 1 _ (w-z) y l + \z\ J 

i i f Vsc-^ 2S ~ 2rn ^\ 

= 1 - j[c ke Ai(s k + si)] 9 a + + -[c ke Ai(s k + s^)]®<7_ + O I y /s + \ x \ J (5 '"^ 

where the notation [Am] stands for a matrix with entries A k g. This yields 

[p x ] M = -i lim AS 12 (A) = -c M Ai(« fc + a t ) + O ( VSc-^ 2S ~ 2m ^) (5.100) 

Poles. It follows from Thm. 5.1 that under the stated conditions /3\ exists and finite for all s = W 
and hence cannot have poles. 

Symmetry. If C = then (for s £ W) r T (A) = — r(— A) (see (5.5)) and the jump matrices M(A) 
then satisfy 

Aft (A) = d z M-\-\)dz , M (A) = l 2r + c^Cc 6 ^® (u +X , + (A) + g^,_ (A)) . (5.101) 



The contours of jump satisfy (also) 7 + = <y_. Then 5 t(A) := 3 T (A) = a^E.(X)a3. This implies 

(0-3IT2CT3 = — CT 2 ) 

B(A) = l+ ai ' g3 + /Wa +... S- 1 (A) = l- aiSg3 + /3l8 ' T2 +... (5.102) 
A A 

~-t(A) = l-^±^ + ... => a 3 S-t ( A)a 3 = l + ^i^^ + ... (5.103) 
A A 

which shows immediately /3i = /?J (as well as cti = — at). 

Formula for the determinant. From Corollary (4.1) we deduce that Dlndet(Id — Ai 2 ?) = 
— 2iTrai and together with eq. (5.55) we have 

D 2 lndet(Id-^|) = ~4Tr(/3^) (5.104) 

from which the formula follows immediately by integration as in the usual Tracy- Widom distribu- 
tion. Q.E.D 

5.3 Noncommutative PXXXIV 

Similarly we can prove 
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Theorem 5.4 Let ai,a 2 be the coefficient matrices in the expansion of the solution T of Problem 


3.2 as in formula (3.15) and define $>{\] 


:=r(A)e e W® CT3 . T/ien 














(5.105) 


Vf.= \ 'v, a 4 


i[tti, ej]®l — 2{&2j e_y}®<7_ — ej®o~+ . 




(5.106) 


r o — - 1 i 

^ +A ( s _i ai ) o 2 J + A 


i[ai,s]-ja'l -s- 
_2iai42[a 2 , s]+2[s, ai]ai-±(ai) 2 l+i[oi,s 




^ 1 07") 

l . ±U ( J 


Denoting by D = Y]j—j d Sj , so that Ds = 


= 1, we have 








D$ = Vb$ , Vb = A 2 l®cr_ - 


2iDai®<7_ — 1®<7-|- = 


-1 " 
A 2 - 2iDai 




(5.108) 


TTie matrices ai,a 2 satisfy the equations (prime denotes action 


o/d; 






a!" = 8i[oi, s]ai 4 


8ai 4 8z[s, 02] 4 6i(a\) 


2 4 4{a' 1; s} 




(5.109) 


a?, = a^ai 








(5.110) 



Remark 5.6 Differentiating (5.109) and using (5.110) one obtains an ODE for the matrix a\ 

a\ v = 6i{a'{, a[} + 8ia[ [s.cn] +8i[o b (7fli] 4 8is[a' 1 ,oi] 4 4{s,a'/} 4 16a[ (5.111) 

Remark 5.7 If r = 1, then we are in the commutative setting and s = s is just a scalar. Then the 
term involving a 2 in (5.109) drops out and we obtain f = d s ) 

a'(' = 8ai 4 6i(a[) 2 4 8s a[ (5.112) 

If we take the equation (5.112) and we differentiate once we obtain, for a' l; the equation (1.19) 
(up to rescaling). For this reason we will call the system (5.109, 5.110) the noncommutative 
Painleve XXXIV equation. 

Corollary 5.2 Denoting by F[ nc ^(s) the Fredholm determinant of the operator Id 4 Ai(»;s) on 

L 2 (R + ,C r ) we have 

d S] F[ nc \s) = l {a l ) 33 (5.113) 
and ai is a solution of the noncommutative PXXXIV equation (5.109, 5.110). In particular 

T)F[ nc) {s) = iTr ai . (5.114) 
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Proof of Cor. 5.2. The Fredholm determinant equals the determinant of Id + K. as explained 
already in the proof of Thm. 5.1. Then the formulas above follow simply from Thm. 4.2 and 
Corollary 4.1. Q.E.D 

Before concluding the paper with the proof of Thm. 5.4 we point out that -in a sense- all the 
relevant information is already contained in Thm. 5.3 because of the matrix Miura relation 



ai (s) Fv ° P = :2 a 1 (s)- i /? 1 (s)° q ^~2i I f3i(t + 5 1 ,...,t + 5 r )dt-if3 1 (s) 



(5.115) 



where s= (S + Si, . . . S + S r ) and {5\ is the noncommutative Hastings-McLeod family of solutions 
(depending on C) in Theorem 5.3. This immediately yields the 



Corollary 5.3 The Fredholm determinant of the matrix Airy convolution kernel Aig satisfies 



det (Id + Aig) = cxp 



Tr (^(t + Si,...,t + S r ) + 2(f - S)ft(t + S x ,...,t + S r )) At 



(5.116) 



where P\(s) is the Hastings- McLeod family of solutions to noncommutative Painleve II as in Thm. 
5.3. 



Proof of Thm. 5.4. Recall that L = 1®L and 
L0-3L- 1 



i 1 
<j + — iXo_ , La 2 L^ 1 = —<J+ + Aer_ , Lo\L~ x = 03 



A 



(5.117) 



L(A) 



dj(ai®a 3 ) 



a 3 T(A) + iXT(\)e j9 a 3 = (v} 2) \ 2 + V^X + vf >) 3(A) 



, ai«(T3 + bi&ai a 2 «>l + &28>0i 
iXL(X) 1 H ! 



(yW\ 3 + vW\ + vj®): 



A 2 

CLi@>0~3 a2®l + 62®Cl 



e,®cr 3 



A 



A 2 



(5.118) 



iX I 1 — zai®cr_ 



(v/ a) A a + vf>A + if>) 



1 

A 2 ' 

<z 2 ® 1 + &2®f3 



A 2 



02®1 + 628173 

A 2 



[ \ a + ~ iXa- 



i 02®1 + &2®CT3 

1 - iai®a- + -^ai»u + H 



(5.119) 



(5.120) 



A 2 (l — iai®o~-)ej®o~- + (iai®<7+ + a2»l + &2®C3) ejgxr- — (1 — iai®CT— )ej®er+ = 

= A 2 vf ] (1 - iai®<7_) + Vj 2) {iai®a + + a 2 ®l + & 2 »cr 3 ) + V". (0) (l - ia ie cr„) (5.121) 
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We thus have 



(<2\ (n\ 

Vj (1 — ia 1 rr :■ = (1 — iai®<j-)ej®<J- => QjUXJ- 

Vj° — i[oi, ej]»l + ([a 2 ,ej] — {&2,ej} — [oi,ej]ai — idjai)»a- — ej»a + 
Looking at the A -1 coefficient one finds the following identity 

-idjai = [ej,a 2 ] + [ai,ej]ai - {e 3l b 2 } 

which allows us to rewrite 



(o). 



£[ai,ey]®l — 2{6 2 ,ej}®cr_ — ej®a+ 



Summing up (5.124) for j = 1, . . . r we also have 

Dai = — 2*62 • 



(5.122) 
(5.123) 

(5.124) 

(5.125) 

(5.126) 



We will need also more information from (5.118) by looking at the coefficients of the negative powers 
of A in particular we will need this for Vd '■= X^=i- To do so it is convenient to multiply (5.118) 
on the left by L _1 . Below we list the results of lengthy but completely straightforward inspections. 
We list the entry of the coefficient of A-' in the form [A J ']fc i £. 



[A- 


-11 
J 1,2 


6 2 = §ai 


[A- 


-2l 
J 1,2 


, 1 / i n 
b 3 = ~2 a i a i ^ 4% 


[A- 


J 1,1 


' 1 ' h 
a 2 — —a x a\ — w 2 a\ = a x a\ 


[A- 


-3l 
J 1,2 


64 = \a[b 2 - 


- \b\ + |aia 2 + \b' z + \b 2 a 2 = -±(ai) 2 + l 2 o! x a 2 - \a![a x - faf 



(5.127) 



A similar and completely straightforward computation (involving longer algebra) yields 



B(X)= 



2 





A(s-6 2 ) 







(5.128) 



63 — ib 2 a\ + i[ai, s] 



-s - 6 2 



iai — 64 - {s, 62} — ib 3 cii + [a 2 , s] + [s, ajai + b 2 (b 2 — af + a 2 ) i[oi,s] + ib 2 a x — o 3 ^^' 
The expression can be simplified using (5.127) to give 



b(x)- 





T +A(s-K) 







i[oi,s] - |< 



2"1 



iai + [s, a 2 ] - §{s,a'i} + [x, ai]ai + ^(a'J 2 + ±a'{' 1 + i[ai,s] + |a" 



(5.130) 
(5.131) 
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One then has to write the zero-curvature equations 



dxV o -DB+[V o ,B]=0 



(5.132) 



which yield 



8i[oi, s]oi + 8a! + 8i[s, a 2 ] + 6i(ai) 2 + 4{a' 1 ,s} 



(5.133) 
(5.134) 



a 2 = ajai 



where the first equation comes from the entry (1, 1) of the coefficient in A~ x of (5.132), while the 
second equation comes from (5.127); all other entries of (5.132) are then automatically zero. One 
can use (5.133) to simplify further the expression for B as given in the statement of the theorem. 
Q.E.D 
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